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Abstract 

The D-dimensional cosmological model on the manifold M = R X Mi x M2 describing the evolution of 2 
Einsteinian factor spaces, M\ and M2, in the presence of multicomponent perfect fluid source is considered. 
The barotropic equation of state for mass-energy densities and the pressures of the components is assumed 
in each space. When the number of the non Ricci-flat factor spaces and the number of the perfect fluid 
components are both equal to 2, the Einstein equations for the model are reduced to the generalized Emden- 
Fowler (second-order ordinary differential) equation, which has been recently investigated by Zaitsev and 
Polyanin within discrete-group analysis. Using the integrable classes of this equation one generates the 
integrable cosmological models. The corresponding metrics are presented. The method is demonstrated for 
the special model with Ricci-flat spaces Mi, M% and the 2-component perfect fluid source. 
PACS numbers: 04.20.J, 04.60.+n, 03.65. Ge 



1 Introduction 

Following the purpose to study the early universe we develop the multidimensional generalization Iql-flR 



|L6|-|21 ,||25|,||26|,||28|1 of the standard Friedman-Robertson- Walker world model. If the extra dimensions of 
the space-time manifold really exist, the unique conceivable site, where they might become dynamically 
important, seems to be possible. This is some early stage of the evolution. Usually within multidimensional 
cosmology (see, for instance, §-0,(l^-[|2^,(2|], j^-|36) and references therein) it is assumed the 

occurrence of the topological partition for the multidimensional space-time on the external 3-dimcnsional 
space and additional so called internal space (or spaces) due to the quantum processes at the beginning 
of this stage. In correspondence with such partition the space-time acquires the topology M = R x M\ x 
. . . x M n , where R is the time axis, one part of the manifolds Mi, . . . ,M n is interpreted as 3-dimensional 
external space and the other part stands for internal spaces. Usually the internal spaces are compact, 
however the models with noncompact internal spaces are also discussed |l4| , [^0| , [^9| |50| . The subsequent 
evolution of the multidimensional Universe is considered as classical admitting the description by means of 
the multidimensional Einstein equations. Achieving the integrability of these equations is the main goal of 
our investigation. As the present world seems to be 4-dimensional, there is the assumption that the internal 
space(s) had contracted to extremely small sizes, which are inaccessible for experiment. This contraction 
accompanied by the expansion of the external space is described by some models (the first model of such 
type has been found in Q) within multidimensional cosmology and is called dynamical compactification. 

We consider a mixture of several perfect fluid components as a source for the multidimensional Einstein 
equations. Such multicomponent systems are usually employed in 4-dimensional cosmology and are quite 
adequate type of matter for description some early epochs in the history of the universe Q] . 

The paper is organized as follows. In section 2 we describe the multidimensional cosmological model and 
obtain the Einstein equations in the form of the Lagrange-Euler equations following from some Lagrangian. 
Here we develop the n-dimensional vector formalism for the integrating of the equations of motion. Conclud- 
ing Section 2 we present the review of the all known integrable models. In Section 3 we suggest the method 
for obtaining the new class of the integrable models on the manifold M — Rx Mi x M^. The method is based 
on the reducing of the Einstein equations to the generalized Emden-Fowler (second-order ordinary differen- 
tial) equation. The method is useful for any 2-component model on the manifold M = Rx Mi x Mi except 
for the cases admitting the integration by more simple way. The total number of the model components is 
equal to the sum of the number of the non Ricci-flat spaces with the number of the perfect fluid components. 



The integrable classes recently derived by Zaitsev and Polyanin of the generalized Emden-Fowler equation 
allow to generate the new integrable cosmological models. Their metrics are presented. In Section 3 the 
method is applied for the models with Ricci-flat spaces Mi, Mi and 2-component perfect fluid. 



2 The model and the equations of motion 

Within n-factor spaces cosmological model Z?-dimensional space-time manifold M is considered as a product 
of the time axis R and n manifolds Mi, . . . , M n , i.e. 

M = R x Mi x . .. x M n , (2.1) 

The product of one part of the manifolds gives the external 3-dimensional space and the remaining part 
stands for so called internal spaces. The internal spaces are supposed to be compact. Further, for sake of 
generality, we admit that dimensions N — dimM^ for i = 1, . . . , n are arbitrary. 
The manifold M is equipped with the metric 

n 

g = - e 27(t)dt®dt + ^2ex-p[2x i (t)]gM, (2.2) 

i=l 

where j(t) is an arbitrary function determining the time t and is the metric on the manifold Mj. We 
suppose that the manifolds Mi, . . . , M„ are the Einstein spaces, i.e. 

Rkidg®} =^41' h,k = l,---,N h i=l,...,n, (2.3) 

where Ai is constant. In the special case, when Mj is a space of constant Riemann curvature i*Q the constant 
Aj reads: A* = K^N - 1) (here Ni > 1). 

Using the assumptions (2.3) we obtain the following non-zero components of the Ricci tensor for the 
metric (2.2) || 

= e" 2 ^ (j2 N ^ l f + 7o - 77'o) (2-4) 
K; = {A* exp[-2^] + [x* + i 4 (7o - 7)] e" 27 } C ( 2 - 5 ) 

where we denoted 

n 

-y = ^2N i x i . (2.6) 

i=l 

Indices m, and rij in (2. 4), (2. 5) for i = 1, . . . ,n run over from (D — fy) to (D — Nj + Ni) 

(D = 1 + £JLi Ni = dimM). 

We consider a source of gravitational field in the form of multicomponent perfect fluid. The energy- 
momentum tensor of such source under the comoving observer condition reads 



■N 



= Y,T^\ (2.7) 



(T%M) = diag (-pl»\t),pM(t)8£,. . . , P M(t)St) > ( 2 - 8 ) 

Furthermore we suppose that for any /x-th component of the perfect fluid the barotropic equation of state 
holds 

p?\t) = (l-hV>)pM(t), A * = l,...,m ) (2.9) 

where = const. It should be noted that each /i-th component admits different barotropic equations of 
state in the different spaces Mi, . . . , M n . From the physical viewpoint this follows from the separation of 
the internal spaces with respect to the external one and with respect to each others. 
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One easily shows that the equation of motion Vm^o^^ = for the fi-th component of the perfect fluid 
described by the tensor (2.8) reads 



1=1 



Using the equations of state (2.9), we obtain from (2.10) the following integrals of motion 



A (m) = p W exp 



= const. 



(2.10) 



(2.11) 



The Einstein equations RfJ — RSff /2 = n 2 T^f (k 2 is the gravitational constant), can be written as 
RM = K 2 [T™ - T5™/(D - 2)]. Further we employ the equations R° - R/2 = k 2 T ° and R£* = k 2 [T^ - 
T5™;/(D - 2)]. Using (2.4)-(2.9), we obtain for them 



n 

- J2 Gi^iP + V = 0, 



(2.12) 



A*e" to< + [x i + x\io - j)]e~^ = -n 2 £ A™ (hf E *g *f 



k 



x exp 



N i Kf ) x i - 2 70 



Li=l 



Here 



Gij = NAj - NiNj 
are the components of the minisupcrspace metric, 

n rh 

-^x'Nie-**' +k 2 J2 aM ex P 

i=l (U=l 



^iV^^-270 



i=l 



(2.13) 
(2.14) 

(2.15) 



The dependence on the densities p^ in (2. 12), (2. 13) has been canceled according to the relations (2.11). 

It is not difficult to verify that after the gauge fixing 7 = F(x 1 , . . . ,x n ) the equations of motion (2.13) 
may be considered as the Lagrangc-Euler equations obtained from the Lagrangian 



(2.16) 



under the zero-energy constraint (2.12). 

Now we introduce n-dimcnsional real vector space R n . By e\, . . . e n we denote the canonical basis in R n 
(ei = (1, 0, . . . , 0) etc.). Hereafter we use the following vectors: 
the vector we need to obtain 

x = x 1 (t)e 1 + ... + x n (t)e n , (2.17) 
the vector induced by the curvature of the space 



2 v— 2 ,* 

1=1 

the vector induced by /x-th component of the perfect fluid 



k 



»=1 



L> - 2 



(2.18) 



(2.19) 
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Let < ., . > be a symmetrical bilinear form defined on R" such that 

<e i ,e j >=G ij . (2.20) 
The form is nongenerated and the inverse matrix to (Gij) has the components 

6 ij 1 

The form < ., . > endows the space R n with the metric, which signature is (—,+,...,+) [p^7|,p8[. By the 
usual way we may introduce the covariant components of vectors. For the vectors Vk and we have 



v w = - 2 §-> «i (fc) =E G *W*) = 2 ( JV *- ,y *)' (2 - 22) 

A i—i 



h (M) _ hl^h , u <f) = ^ Gyti^) = A^- (2.23) 



The values of < Vk, >, < Vk, Un > and < Un,u u > are presented in Table 1. 



< .,. > 


V 3 




Vi 


4(fe-l) 


-2h^ 


Ufj, 


-2hf 





TABLE I. Values of the bilinear form < ., . > for the vectors Vi and it M , induced by curvature of the space 
Mi and /i-th component of the perfect fluid correspondingly. 

A vector y 6 R n is called time-like, space-like or isotropic, if < y, y > has negative, positive or null 
values correspondingly. Vectors y and z are called orthogonal if < y, z >— 0. It should be noted that the 
curvature induced vector Vi is always time-like, while the perfect fluid induced vector w M admits any value 
of < 



> (see Table 1). 



Using the notation < ., . > and the vectors (2.17)-(2.19), we may write the zero-energy constraint (2.12) 
and the Lagrangian (2.16) in the form 



1 

E = - < x,x > +V = 0, 



L = c 70 " 7 | < x,x > -V 



(2.24) 
(2.25) 



where 



y = c 2 (7-7o) 



2 ^ 

i=l 



+ K 



E^ ( 



(2.26) 



It is obviously from (2.26) that the term induced in the potential by the non-Ricci flat space Mi is similar 
to the term induced by /i-component of the perfect fluid. Due to this fact the non-zero curvature of the 
manifold Mi may be also called a component and now we use the notion of the component in such new sense. 
Further we employ the so called harmonic time gauge, which implies 



7(t) = 7o = E N « 



(2.27) 



4 



^From the mathematical viewpoint the problem consist in integrability of the system with n > 2 degrees 
of freedom, described by the Lagrangian of the form 

L = - < ±,x > -^oMe^"^, (2.28) 
n=x 

where x, G R n . In (2.28) m denotes the total number of the components including the curvatures and the 
perfect fluid components. It should be noted that the kinetic term < x, x > is not positively definite bilinear 
form as it usually takes place in classical mechanics. Due to the pseudo-Euclidean signature (— , +, +) of 
the form < ., . > such systems may be called pseudo-Euclidean Toda-like systems as the potential like that 
given in (2.28) defines well known in classical mechanics Toda lattices |54|. In the papers p[,p9|,p0[ the 
following classes of the integrable pseudo-Euclidean Toda-like systems have been found 

1. m — 0. This case corresponds to the vacuum multidimensional cosmological model on the manifold 
M = Rx Mi x ... x M n with all Ricci-flat spaces M*. The corresponding metric is a multidimensional 
generalization of the well-known Kasner solution jl9j . 

2. m=l, the vector b\ is arbitrary. The metrics for this 1-component case were obtained in po|| . This 
integrable class may be enlarged by the addition of the new components inducing the vectors collinear 
to the vector b\. 

3. m > 2, n = 2, b^ = b + C^bo, where b is an arbitrary vector and bo is an arbitrary isotropic vector, 
Cp=const. This class was integrated in [^0| only under the zero energy constraint. 

4. m > 2, the vectors b\, . . . ,b m are linear independent and satisfy the conditions 
< b^^bv >= for fi =/= v. This integrable class may be enlarged by the addition of the new com- 
ponents inducing the vectors collinear to one from the orthogonal set b\,...,b m . The corresponding 
cosmological models are studied in || . 

5. m > 2, the vectors b%, . . . ,b m are space- like and may be interpreted as a set of admissible roots || 
of a simple complex Lie's algebra G. In this case the pseudo-Euclidean Toda-like system is trivially 
reducible to the Toda lattice associated with the Lie algebra G j34|. For G = A2 = sl(3, C) the metric 
of the corresponding cosmological model was explicitly written in Q . 

In the present paper we consider only 2-component (m = 2) pseudo-Euclidean Toda-like systems with 2 
degrees of freedom (n = 2) under the zero energy constraint. The corresponding multidimensional cosmo- 
logical models are 2-factor spaces, i.e. 

M = Rx Mi x M 2 (2.29) 

and admit the following combinations of the components: curvature of Mi and curvature of M2 (vacuum 
models); curvature of Mi or Mi and 1-component perfect fluid; 2-component perfect fluid in Ricci-flat 
spaces Mi and M 2 . In our recent paper we have integrated the vacuum model of the type (2.29) with 
2 curvatures for the dimensions (dimMi,dimM 2 )=(6,3),(8,2),(5,5). Now we develop more general procedure 
useful for any combination of the 2 components. 



3 Reducing to the generalized Emden- Fowler equation 

Let us consider the equations of motion following from the Lagrangian (2.28) with n = m = 2 under the 
zero energy constraint. If the vectors b\ and 62 satisfy one of the following conditions 

1. b\ and 62 are linearly dependent, 

2. < bi, b 2 >= 0, i.e. b\ and &2 are orthogonal , 

3. < b\ — 62, b\ — 62 >= 0, i.e vector b\ — bi is isotropic, 
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the equations of motion are easily integrable and the corresponding exact solutions have been obtained in 
the papers ||,|2(J. Now we aim to develop the integration procedure just for all remaining cases. Then, 
further we suppose that the vectors b\ and b 2 do not satisfy any condition from 1-3. 
Let us introduce in R 2 an orthogonal basis forming by the following two vectors 



h = ("22 - "12)61 + ("11 - u 12 )b 2 , h = b2 - h, 



where we denoted 



Upv =< bfj,,b u >, j«,i/=l,2. 
According to the admission accepted f 2 is not isotropic vector, i.e. 

< h, j*2 >= "11 + "22 - 2lti2 7^ 0. 



(3.1 
(3.2 

(3.3 



One may easily check that u\ 2 ~ "n"22 > for any vectors b\, b 2 £ R 2 and u 2 2 — unu 22 — if and only if 
b\ and b 2 are linearly dependent. Then in the case under consideration 



< fl,fl >= -("12 - "ll"22)("ll + "22 - 2Mi 2 ) ^ 0, 
< fljl > I < h,f2 >= -{u\ 2 - "ii" 22 ) < 0, 

i.e. one from the orthogonal vectors /1 and f 2 is space-like and the other is time-like. 
The vector x(t) we have to find decomposes as follows 

< x, ft > < xj 2 > 

~h + ——? — T^rte- 



< h,h >' 



< h,h >' 



For the new configuration variables 



z(t) = 
1 



< x, h > 



In 



1.(2) 



y(t) = 



< h,h > 

< fufi > 



<xji> 



the Lagrangian (2.28) and the corresponding zero-energy constraint look as follows 

L = 2p(z 2 -y 2 ) -V(z,y), 
E = 2f3(i 2 -y 2 ) +V(z,y) = 0, 

where the potential V{z, y) has the form 



V(z,y) = V e 2al3y (sgn 



sgn 



,(2) 



J2f3 2 z 



In formulas (3.9)-(3.11) the following constants are used 

a = V u i2 - "n"22, P = ("11 + "22 - 2mi 2 ) _1 , 

01 = -("11 - Uu)(3, 02 = P\ + 1 = ("22 - "l2)/3, 



(3.4 
(3.5 



(3.6 

(3.7 

(3 



(3.9 
(3.10 



(3.11 



(3.12 
(3.13 
(3.14 



It should be mentioned that using of a basis in the form (3.1) provides the factorization of the potential 
(3.11) with respect to the coordinates of the vector x(t) (the additional linear transformation (3. 7), (3. 8 
does not matter in this situation). Such factorization of the potential is essential under the developing of 
the following procedure proposed in Using the equation of motion following from the Lagrangian (3.9) 



-J-l/ e 2a/3y f/3isgn L«l e 2/3lZ + /? 2 sgn 
lo V L J 



2/3 
! -V(z,y), 



,(2) 



3 2/3 2 z 



(3.15 
(3.16 
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the zero-energy condition (3.10) written in the form 



and the relation 



1 V(z,y) 
W (y/i) 2 - 

_ 

dz 2 



1 V(z,y) 
2/3 (dy/dz) 2 - 



'dz 



we obtain the following second-order ordinary differential equation 

2 



fy 

dz 2 



1 



where 



sgn 



A + /3 2 + 



dy 
dz 



a/3 



(3.17) 



(3.18) 



(3.19) 



(3.20) 



We notice that due to the factorization of the potential the right side of the equation (3.19) does not contain 
y, so, in fact, the equation is the first-order one with respect to dy/dz. 

This procedure is valid for the solutions such that z ^ 0. Under the zero energy constraint the solutions 
of (3. 15), (3. 16) with z = gives the following vector x(t) 



x(t) =p\n\t - t \ + q, 



(3.21) 



where the constant vectors p,q 6 R 2 are such that 



P = —oh, 
a z 

„<?,&!> _ 



02 



a«a 2 /3 



> 0, e 



<<j,fc 2 > _ 



01 



a( 2 )a 2 /3 



> 0. 



(3.22) 
(3.23) 



We note that the exceptional solution (3.21) exists only if the inequalities in (3.23) are satisfied. It should 
be mentioned, that the set of the equations (3. 10), (3. 15), (3. 16) does not admit static solutions z = y = 
due to the condition (3.3). The solutions with z = ±y are also impossible, so using the relation (3.17) we 
do not lose any solutions of the set (3. 10), (3. 15), (3. 16) except, possibly, the solution (3.21). 

Let us suppose that one is able to obtain the general solution of the equation (3.19) in the parametrical 
form z = z(r), y — j/(t), where r is a parameter. Then using (3.6)-(3.8) we obtain the vector x as the 
function of the parameter r 



x{r) ^'^^(-(3 2 b 1 +l3 1 b 2 ) + 2(3 



z(t) — In - 



7,(1) 



(6 2 -6l). 



(3.24) 



We recall that coordinates of the vector x(t) in the canonical basis are the logarithms of the scale factors for 
the spaces Mi, M 2 . The relation between the harmonic time t and the parameter t may be always derived 
by integration of the zero-energy constraint written in the form of the separable equation 



dt 2 = 2/3 



fdy 



dzV 
dr) 



V(z(r),y(r)) 



dr 2 . 



(3.25) 



Thus the problem of the integrability by quadrature of the pseudo-Euclidean Toda-like systems with 2 degrees 
of freedom under the zero-energy constraint is reduced to the integrability of the equation (3.19). 

For dy/dz the equation (3.19) represents the first-order nonlinear ordinary differential equation. Its right 
side is third-order polynom (with the coefficients depending on z) with respect to the dy/dz. An equation 
of such type is called Abel's equation (see, for instance |27]] , |37[] ) . There are no methods to integrate 
arbitrary Abel's equation, however the equation (3.19) may be integrated for some values of the parameters 
a/3 and /3i + f3 2 - First of all let us notice that the equation (3.19) has the partial integrals y ± z =const, 
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which make the relation (3.17) singular and as was already mentioned are not partial integrals of the set 
(3. 10), (3. 15), (3. 16). Existence of this partial solution of the Abel equation (3.19) allows to find the following 
nontrivial transformation 



,2z 



X dY 
Y dX' 





Y 




z + In 




+ lnC 




X 





5 = ±1, C>0, 



which reduces the Abel equation (3.19) to the so called generalized Emden- Fowler equation 

d 2 Y ( dY 

dA 2 ~ \dX 

where the constant parameters n, m and I read 



- (Pi + pa - 2da/3 - 3) = ■ . 

2 Un + U22 - 2U\2 



1 



m 



I 



(0! +{3 2 - 2Sa/3 + 3) 



-un - 2u 22 + 3ui2 + 8\Ju\ 2 - unu 2 2 
un + u 2 2 - 2ui2 



- (Pi + P2 + 2da/^ - 3) = ■ . 

2 un + U22 - 2ui2 



For our models the parameters in the generalized Emden- Fawler equation are not independent, 
from (3. 29), (3.30) that 

n + m = —3. 

In the special case I — the equation (3.28) is known as the Emden-Fowler equation. 

If the parameters I and m given by (3. 31), (3. 30) are such that 1 = 0, m ^ 1 there exists 
transformation 



1 + ec 2z = 



y 



2 X dY 
' m - 1 Y dX' 

' ■ In Y 2 + C 



m — 1 



, S = ±l, C E R, 



which reduces the Abel equation (3.19) to the following integrable Emden-Fowler equation 

d 2 Y 



dX 2 



m+3 

y—. 



(3.26) 
(3.27) 

(3.28) 

(3.29) 
(3.30) 
(3.31) 
It follows 
(3.32) 

one more 

(3.33) 
(3.34) 

(3.35) 



There are no methods for integrating of the generalized Emden-Fowler equation with arbitrary independent 
parameters n, m and However, the discrete-group methods recently developed by Zaitsev and Polyanin 
[p7| allows to integrate by quadrature 3 two-parametrical classes, 11 one-parametrical classes and about 90 
separated points in the parametrical space (n, m, I) of the generalized Emden-Fowler equation. For instance, 
the two-parametrical integrable classes arise when m and I are arbitrary and n = or when n and I are 
arbitrary and m — 0. The one-parametrical class with / = and n+m = —3 is also integrable by quadrature. 

Let us suppose that the two components of the 2-factor spaces cosmological model under consideration 
induce such vectors b\ and 62 that the corresponding to the model generalized Emden-Fowler equation 
(3.28) with the parameters defined by (3.2),(3.29)-(3.31) is integrable in the parametrical form X = A(r), 
Y = Y(t), where r is a parameter. Then, using the parameter r as the new time coordinate we obtain by 
the formulas (3.26), (3.27), (3.24), (3. 25) the following final result for the metric (2.2) 



-/ 2 (r)[ax(r)] 2 ^ [a 2 (r)} 2 ^dr ® dr + [ ai (r)] 2 g^ + [a 2 (r)] 2 5 (2) , 



(3.36) 
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where we denoted 



f(r) 



l2\/3\ 



(Jn+l 



ai(r) = e 



= « W = e 7* 




[X'(r)} 2 
X(r)Y(r)Y'(rY 



2/3 



By 7* for i = 1, 2 we denoted the following constants 
' LnC 



y( T ) 



7* = 2/? ■ 



n + 1 



(n-l + A)b l (1) -(n-l + 2)b\ 2) 



X(r) 



In • 



(2+n)6J 1) -(l+n)6J 2) ^| 7+^ 



(3.37) 
(3.38) 



,.(2) 



,.(1) 



0(2) - 0(1) 



(3.39) 



We recall that K * are coordinates of the vector 6 M in the canonical basis. In the special case I — one 
may also use by the similar manner the transformation (3. 33), (3. 34) and the result of integration of the 
equation (3.35) to write the metric. This transformation was used in JTo[ | for integrating of the models with 
two curvatures. 

Thus the method described allows to integrate the cosmological models if the corresponding generalized 
Emden-Fowler equation is integrable. Note that if the model with some vectors b\ and 62 is integrable by 
such manner then any model with the vectors ab\ and abi [a is an arbitrary non-zero constant) is also 
integrable as the parameters n, m and I do not change under such transformation of the vectors. Taking into 
account the classes 1-4 (the class 5 does not arise for n = 2) and the additional to them class, which may 
be integrated by the method described, we obtain the quite large variety of the integrable 2-factor spaces 
cosmological models with 2 components. 



4 Examples of the integrable models 

Now we apply the method proposed in Section 3 to the cosmological models on the manifold (2.29) with 
both Ricci-flat spaces M\,Mi and the 2-component perfect fluid source. Let us represent such model by 
Table 2 



manifold / source 


external space M^ 1 


internal space M^ 2 


l-st component of 
the perfect fluid 






2-nd component of 
the perfect fluid 




/4 2) 



TABLE 2. Representation of the model on the manifold M — Rx Mi x M% with Ricci-flat spaces M\, Mi 
for the 2-component perfect fluid. 

We recall that Nj =dimMj and are the constant parameters in the barotropic equation of state 
(2.9). The model is entirely defined by these 6 parameters. One easily shows |l2j that the dominant energy 
condition applied to the stress-energy tensor (2.8) implies < < 2. Usually rational values of the 
parameter are employed in cosmology, for instance, — (iVj — l)/iV, - radiation, = 1 - dust, 
= - Zeldovich (stiff) matter, = 2- false vacuum (A-term), = (D — l)/D - superradiation etc. 
On the other hand the most known cases, when the generalized Emden-Fowler equation (3.28) is integrable, 
arise for the rational parameters n, m and I. So if one demands the rationality of the parameters n, m and / 
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in the equation (3.28) corresponding to the model under the condition of the rationality for the parameters 
, then due to the following relation 



a 2 = u[ 2 - U11U22 = 



N X N 2 



N 1 +N 2 



,(2) 



(4.1) 



the dimensions Ni, N 2 with integer value of the expression R = \J N\N 2 (Ni + N 2 — 1) are singled out. For in- 
stance, the expression R is integer for the following dimensions: (Ni, N 2 ) = (3, 6), (2, 8), (5, 5), (7, 8), (3, 25), (Ni, 1). 
From the physical viewpoint the following cases may be of interest: (2, 1), (3, 1), (3, 6), (3, 25). 

Let us consider the models of the type represented in Table 2 leading to the generalized Emden-Fowler 
equation (3.28) with 

I = 0. (4.2) 

Due to the relation (3.32) arising for our models the equation is integrable for arbitrary parameter m and its 
exact solution has been written by Zaitsev and Polyanin |27j] . It is worth to mention that other 2 integrable 
classes of the generalized Emden-Fowler equation (3.28), arising when n — or m = 0, describe the same 
cosmological models. It follows easily from (3. 29), (3. 31) that if the model is such that I = for 5 = 1 (or 
5 = —1) then n = for 5 = — 1 (correspondingly, 6 = 1). It is easy to see also from (3. 30), (3. 31) that 
the condition I = transforms to the condition m = under the inverse numbering of the components. 
Thus, from 3 integrable classes, arising for n = 0, m = and I = 0, correspondingly, of the generalized 
Emden-Fowler equation (3.28) for our models it is enough to study any one from them, let it be the class 
with I = 0. In this case the equation (3.28) has the form 



d 2 Y 
dX 2 



X- m - 3 Y m , m = -[3i - fa 



By the following transformation 
it reduces to the equation 



which is easily integrable. Then the general solution of the equation (4.3) has the form 



Y=±* X=± ' 



■F{ry 



F(rY 



where 



F(r) - / 



771+1 



m+1 



•Ci 



-1/2 



dr + C 2 , m ^ — 1, 



y [2 ln|x| +dr 1/2 dT + C 2 , m = -l. 



(4.3) 
(4.4) 
(4.5) 

(4.6) 

(4.7) 
(4.8) 



We suppose that the both components of the perfect fluid have the positive mass-energy densities given by 
(2.11). It means = k 2 A^ > for /j, = 1, 2, so e = sgn [aWe^ 2 )] = 1. Then taking into account the 
formula (3.26), we must consider the general solution (4.6) on such interval of the variable r where 



F(r) 

Y(t)X'(t) tF'(t) 



- 1 > 0. 



(4.9) 



Finally using the results of Section 3 we obtain the following exact solution for the cosmological model 
represented by Table 2 in the special case / = 0: 
the metric is given by the formula (3.36), where 



/ 2 (r) 



2\(3\C- 25a0 



F'(r) 



G(t)t 2 



(4.10) 
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the scale factors of the spaces Mi , M 2 have the form 



«i(r) 



A( 2 ) 



[G(7 



2 T 2l/ 3 i«j2)-/ 3 2«(l) 



5/a 



the mass-energy densities of the components read 



P (1) W = 



^ (2) (r) 



[ai(r)] 2 ^[a 2 (r)] 
[A«f 2 [AW]" 



2W 2 



{[G(r)p 
{[G(r)} 



2— 2«n 



2 T 2j Q 0j 



5/a 



5/q 



[ai(r)] 2 ^[a 2 (r)] 2 ^ 
The functions F(r) and G(r) are defined in (4.7)-(4.9); components 

are given in (2.23); the parameters a,/?, /?i,/?2 are defined in (3. 12), (3. 13); the values u^ v =< u 
fi,u = 1,2 may be calculated from the parameters Ni, fti by the formula given in Table 1. 
The following relation is valid for the densities (4.12)-(4.13) 



(4.11) 

(4.12) 
(4.13) 



2_2j a 0" 

of the vectors in the canonical basis 

, > for 



(4.14) 



We recall the possible existence of the special solution (3.21). 

Concluding the paper, we mention in Table 3 some interesting from our viewpoint special models with 
I = for the dimensions Ni = 3 and _/V 2 = 6. One easily shows that for these dimensions and given in 
Table 3 values of the parameters the parameter I given by (3.31) is equal to zero, so the special model 
is described by the exact solution (4.10)-(4.13). It follows from (2. 22), (2. 23) that for these dimensions the 2 



perfect fluid components with the parameters 1-st: h\ x> = 4/3. /A = 2, 2-nd: h\ = 2. h 2 = 5/3 induce 
the vectors, which coincide with the vectors induced by the curvatures of Mf and Mf , correspondingly. 
Then the adding of such 2 components to the integrable vacuum model on R x Mf x Mf with 2 curvatures 
(see investigation of this model in ficfl ) provides with the integrable model for the 2 non Ricci-flat spaces 
and the 2-component perfect fluid. 



manifold/source 


external space Mf 


internal space Mf 


l-st component of 
the perfect fluid 


radiation 

h {1) - * 
'h — 3 


radiation 

h {1) - £ 

n 2 - fi 


2-nd component of 
the perfect fluid 


radiation 

h {2) - 1 
'h ~ 3 


Zeldovich matter 

/4 2) = 


l-st component of 
the perfect fluid 


dust 
= 1 


radiation 
h {1) - ^ 


2-nd component of 
the perfect fluid 


radiation 

h {2) - 1 
"1 — 3 


dust 

/4 2) = 1 


l-st component of 
the perfect fluid 


radiation 

h {1) - 1 
"1 — 3 


radiation 
h {1) - ^ 

n 2 ~ fi 


2-nd component of 
the perfect fluid 


dust 

/4 2) = 1 


radiation 

"2 — 6 


l-st component of 
the perfect fluid 


radiation 

h {1) - 2 
n i — 3 


radiation 
h {1) - ^ 

"2 — 6 


2-nd component of 
the perfect fluid 


false vacuum 
h{ 2) =2 


false vacuum 

/4 2) = 2 


l-st component of 
the perfect fluid 


false vacuum 
fcW - 2 


false vacuum 
tf> = 2 


2-nd component of 
the perfect fluid 


Zeldovich matter 
h^=0 


false vacuum 

/4 2) = 2 
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TABLE 3. Examples of the integrable models for dimensions Ni = 3 and N 2 — 6. The corresponding 
exact solutions are given by the formulas (4.10)-(4.13). 
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